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Abstract 

We derive equation describing distribution of energy losses of the 
particle propagating in fractal medium with quenched and dynamic 
heterogeneities. We show that in the case of the medium with frac- 
tal dimension 2 < D < 3 the losses of energy are described by the 
Mittag-Lemer renewal process. The average energy loss of the parti- 
cle experiences anomalous drift A ~ x a with power-law dependence 
on the distance x from the surface and exponent a = D — 2. 
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1 Introduction 



In this paper we study the problem of energy losses of the particle propagating 
in disordered medium. We consider both a model of discrete energy loss, 
when the particle collides with randomly placed scattering centers in the 
medium, so that during intervals between such scattering events the particle 
has ballistic trajectory and the energy of the particle does not change, and its 
continuum limit, corresponding to the description of the energy loss process 
in terms of kinetic equation. 

There exists a close analogy between this problem and the diffusion of a 
particle through the medium: the energy loss is proportional to the number 
of scattering events n (x) at the distance x from the surface of the medium, 
and x plays the role of the time t in the diffusion problem. The propagation 
of a particle in the regular homogeneous medium, with the number of scat- 
tering events (n) ~ x, is analogous to the standard Brownian motion with 
the average number of events (n) ~ t. Anomalous power-law coordinate 
dependence (n) ~ x a is related to the presence of power-law spatial correla- 
tions in the density distribution of scattering centers in the random medium. 
Such correlations are typical both for fractal systems with quenched disorder 
and for systems with dynamic heterogeneities at the critical point of a phase 
transition. 

In diffusion problems anomalous time dependence of the number of events 
is often observed near glass transition [UI2], for sub-diffusion-limited reactions [3], 
in turbulent regime [1] and for some transport problems in disordered systems. [5], 
i 

Our consideration of the discrete energy loss model with power law cor- 
relations is done at the level analogous to the Continuous Random Walk 
(CTRW) j7]. The main distinction of the stochastic process describing the 
energy loss and the CTRW one is that the energy loss problem corresponds to 
a sum of random positive quantities - energy losses at scattering events. Of 
substantial interest is a continuous limit of the discrete energy loss model cor- 
responding to some effective kinetic description. The highly nontrivial and 
interesting transition from the discrete to kinetic description in the usual 
CTRW case was discussed in [8]. 

In section [2] we study the energy losses in weakly heterogeneous medium. 
We demonstrate that average energy losses in such medium can anomalously 
grow with the distance x from the surface, (A(x)) ~ x 1///3 , with super- 
diffusional exponent (3~ l > 1. In section[3]we study energy losses in strongly 
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heterogeneous medium. We first formulate simple model of anomalous scat- 
tering, describing the loss of energy by the particle during its propagation in 
random medium (section [3. ip . In section [3721 this process is generalized to the 
case of general distribution of energy losses during individual collisions using 
formalism of CTRW. We derive equation for the distribution function of en- 
ergy losses and show that the energy of the particle experiences anomalous 
drift with average energy losses growing with the distance x as (A(x)) ~ x a , 
with sufe-diffusional exponent a < 1. In section [3731 we demonstrate, that the 
loss of the energy because of scattering on fractal structure is described by 
the Mittag-Lemer renewal process [9], which may be considered as fractional 
generalization of the well known Poisson renewal process. Main results of 
this work are summarized in section [H 



2 Weakly heterogeneous medium 

The loss of the energy by the particle propagating in homogeneous medium 
is described, in the continuous limit, by the Landau kinetic equation for the 
distribution function of energy losses at point x\10\ 

df (A,x) _ 1 
dx a 



POO 

/ dew(e)[f{A-e,x)-f{A,x)} (1) 
Jo 



where w (e) is the probability distribution of energy loss at scattering event 
and a -1 is linear density of scattering centers separated by a distance a. 

In the random medium the probability distribution of energy loss becomes 
a random function of the coordinate x, w (e) — > w(e\x). In principle one 
could imagine arbitrary variations of the form of w (e\x) from one point to 
another that could be correlated over spatial domains with the size controlled 
by the corresponding correlation length £. 

In the simplest case of finite correlation length £ the basic kinetic equation 
for the averaged distribution function / (A,x) on distances x ^> £ takes the 
form 

df (A,x) 1 
dx 

with 



POO 

- dew{e) [f (A — e,x) — f (A, x)) (2) 
'Jo 



w(s) = (w(e\x)}. (3) 
The averaging in Eq. ([3]) is over the randomness of the distribution function 



W[£\X) 
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As an illustration, let us consider the distribution w(e\a) used in Eq. ([T]), 
where a denotes the set of parameters characterizing this distribution. The 
simplest way of introducing spatial randomness of the distribution w(e\a) is 
to consider the parameters a as random fields depending on spatial coordi- 
nates, w(e\a) — > w(e\a(x)). At distances larger than the correlation length 
£ the fields w(e\a(x)) are effectively uncorrelated and the distribution w(e) 
is simply given by 

w(e) = J dap(a)w(e\a) (4) 

where p(a) is the probability distribution describing the local fluctuations 
of cr(x). The properties of the resulting distribution w(e) can substantially 
differ from those of the initial distribution w(e\a). Of special importance is 
the influence of randomness in a on the asymptotic behavior of w(e) at large 
e. In particular, the power-like decay of p(a) at large a induces the power-like 
decay of w(e) even for well-localized exponentially decaying w(e\a). 

The solution of Eq. ([2]) is most conveniently obtained by solving for the 
Laplace-transformed distribution function 

POO 

f(p, x) = / dAe~ pA f(A, x) = exp [-co(p)x] (5) 
Jo 

where 

1 f°° 

u(p) = - / dew(e) (l - e~ pe ) (6) 
a Jo 

Calculating the inverse Laplace transform of Eq. (jSJ) we find that the 
distribution function / (A,x) can have different form depending on the large 
energy asymptotes of the distribution w (e): 

a) If this function w (e) decays at large e faster than 1/e 3 , the distribution 
function / (A,x) has Gaussian form with the center at A = Aix and the 
width y/A~2X, where are corresponding moments of the distribution w (e): 

POO 

A k = / dee k w(e) (7) 
Jo 

b) If the function w (e) decays as l/e 1+l3 with 1 < < 2 the distribution 
function / (A,x) is still centered at A = A\x but is not Gaussian, and its 
width grows as power law x 1 ^ of x. In both cases a) and b) the distribution 
/ (A,x) becomes sharper with the rise of x. 
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c) If the function w (e) decays slower than 1/e 2 then the first moment 
Ai diverges. In the case of power distribution w (s) ~ e^ 1 ^ 13 with exponent 
< f3 < 1 we find from Eq. (0) u (p) — Cp 13 with certain constant C ~ A 13 /a. 
Calculating the inverse Laplace transformation of Eq. (jSJ) we get 

where Wp (z) is the Wright type function (see Appendix |B]). The center 
of this distribution and its width are on the same order of value and grow 
as the power law (A (x)) ~ x 1 ^ of the distance x from the surface. Such 
dependence is familiar for Levy flight processes characterized by the exponent 
> l.PU 

3 Strongly heterogeneous medium 

In this section we consider energy losses in strongly heterogeneous medium 
with infinite correlation radius £ (or on distances from the surface small with 
respect to £). Such situation takes place in the case of fractal medium or for 
the system at critical point of phase transition. 

3.1 Simple model with constant losses 

Let us first consider the simple one-dimensional microscopic model for energy 
loss in heterogeneous medium in which the loss takes place at through a 
sequence of discrete events in which the incident particle looses an equal 
amount of energy A x . This model describes the eikonal energy loss of a 
high energy particle loosing energy in small portions so that one can, in 
the first approximation, consider its trajectory to be a straight line. The 
model generalizes the simplest model of discrete random walk with the fixed 
elementary step to the case of the positive variable. The assumption of 
constant energy loss will be relaxed in the next paragraph in which we will 
consider a more general model. More precisely, in the simple model we 
consider a particle entering the medium at the point x = and loosing 
energy at random points Xi, ■ ■ ■ ,x n lying on the trajectory composed by the 
intervals of length l 1: • ■ ■ ,l n , see Fig. [TJ 

The cumulative energy loss at some point x for the given event is, evi- 
dently, equal A (x) = Ain (x), where n(x) is the number of scattering events 
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Figure 1: Particle propagating in random medium looses the energy Ai dur- 
ing scattering events which are separated by random intervals I along the 
trajectory of the particle. 



that the particle experienced along its trajectory from the initial point to 
the point under consideration. The probabilistic description of the process 
of random energy loss is thus fully specified by the probabilistic properties 
of n(x). Let consider the case in which these properties are fully described 
by the probability density if) (I) of the spatial distance I = x — x' between 
the points x and x' < x at which two subsequent scattering events took 
place and the distribution (p(h) characterizes the probability density of the 
position l\ of the first scattering event, see Fig. [TJ In what follows we will 
simply assume that ip(h) = if>(h)- The function if) (I) determines the survival 
probability ^ (r) that the particle experiences no collisions at the distance 
r = x — x' after the collision at point x'\ 

/oo 
if)(l)dl (9) 

To develop a quantitative description of the energy loss process it is con- 
venient to introduce the probability density if) n (x — x') for the n-th collision 
to take place at the distance x — x' from the collision at point x'. Then 
the probability density f n {x) for the number n of scattering events along the 
trajectory leading to the point x reads 

f n (x)= [ dx'^(x-x')lf) n (x') (10) 

Jo 

In turn, the probability density if) n (x) is determined by recurrence relation 

/•X 

if) n (x) = / dx'if) (x - x) if) n -i {x') (11) 



o 



with if>i (x) — if) (x). 

The equations fllOflip are easily solved by using the Laplace transform 



fn(q)= / f n (x)e-**dx 
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We have 

fn (q) = r (q) (12) 

q 

where ip (q) is the Laplace transform of the probability density tp(l) . 

Of special interest is the case of the medium characterized by power-law 
scale-invariant fluctuations of the positions of the scattering centers. An ex- 
ample of this situation is provided by the fractal medium. In the considered 
one- dimensional case the relevant correlations are determined by projecting 
the full three-dimensional correlation pattern on the axis of particle propa- 
gation. In the scale-invariant case the function should read 

^ W = { cJ/l^ a! !>a' ^ 

where, as shown in Appendix |AJ exponent a is related to the fractal dimen- 
sion D of the scattering medium as: 

a = D-2. (14) 

Notice, that there are both upper and low boundaries of the fractal dimen- 
sion D for the problem of anomalous scattering. Upper fractal dimension 
D = 3 describes homogeneous medium and at lower critical fractional di- 
mension D = 2 the scattering cluster is too sparse to effectively scatter the 
propagating particle. 

At distances much larger than the ultraviolet cutoff a the Laplace trans- 
form of the function in Eq. ffTB"]) reads 

POO 

4>(q)= ip (x) e~ qx dx ~ 1 - (aq) a , aq < 1 (15) 
Jo 

In the case a < 1 the average scattering length I diverges leading to anoma- 
lous losses of the particle energy. Calculating the inverse Laplace transform 
of Eq. (fl2|) . we get at large x ^> a the distribution function of energy losses: 

f n (x) = (a/x) a W a [(a/x) a n], (16) 

where W a (z) is the Wright type function, defined in Eq. ( HP]) of Appendix [Bl 
Using this distribution, we can calculate the average energy loss at the dis- 
tance x from the surface, which is monotonically increasing function of x: 

(A(x))=A 1 n(x)= ^ 1 (-)", 0<«<1 (17) 
1 1 + a) \a/ 
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At a = 1 the average scattering length I is finite and in average the losses of 
particle energy (A(x)) = Aix/l grow linearly with the distance x, similar to 
the case of homogeneous system [10] . 

3.2 Continuous Time Random Walk (CTRW) 

In this section we consider the case, when energy loss e at each collision with 
scattering center is random and is described by the distribution function 
w (e). By analogy with CTRW [7] the process of energy losses can be con- 
sidered as subordinated to the anomalous scattering process shown in Fig. [TJ 
Let us stress once again that the model under consideration differs from the 
usual CTRW in considering the random walk of the positive variable e. The 
distribution function of energy losses can be presented as the convolution of 
the distribution function of the number of collisions /„ (x) studied in previ- 
ous section and the distribution function of energy loss w n (A) for the given 
number of collisions n: 

/ (A,*) = J>„ (A) /„ (x) . (18) 

n 

Due to the additive nature of the energy loss the distribution function w n (A) 
is determined by the inverse Laplace transform of the power w n (p) of the 
Laplace transform w (p) of the distribution w (e). Explicit expression ffT8]) 
for the distribution function / (A,x) is not very convenient to use, so below 
we derive the kinetic equation for it. 

Simple probabilistic consideration of the losses of the energy by the par- 
ticle propagating in random medium characterized by the probability distri- 
bution w (e) shows that this random process can be described by equation 

px pOO 

f (A,x) = 5 (A) $ (x) + / dx'ip (x - x') dew (e) / (A — e,x') (19) 

Jo Jo 

where $ (x) is the survival probability at the distance x from the surface, 
Eq. (jUJ). To solve Eq. f fl9|) it is convenient to introduce the double Laplace 
transform of / (A,x) over both variables A and x 

poo poo 

/(p,g)= / dAe- pA / dxe-**f(A,x) 
Jo Jo 

The corresponding equation for / (p, q) reads: 

/ (P, q) = — ^ +^{q)w (p) / (p, q) (20) 
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Introducing the function 



H,) = M-, (2D 
1 - -0 (g) 



we can rewrite the above equation (12T|) in the form 

/ (P, 9) = - + 9 (?) [w (p) - 1] / (p, ?) (22) 
Q 

which is equivalent to 

f(A,x) = 5(A) + dx'g(x-x')x (23) 
Jo 



/■oo 

/ dew(e)[f(A-e,x')-f(A,x')} 
Jo 



Eq. ( 123!) can be considered as generalization of Eq. ([I]) for homogeneous 
medium corresponding to constant density of scattering centers, g(r) = 1/a. 
The function g (r) is found by inverse Laplace transform of the function g (q), 
Eq. (1211) . Using Eq. (TI2"]) <? (g) can be related to the Laplace transform of the 
average number of scattering events: 

**(?)=£ nf n (q) = ^-, (24) 

z — 'n g 

We conclude, that g (r) has the meaning of average density of scattering 
events along the direction e x of particle propagation at the distance r from 
the scattering event: 

g (r) = dn (r) jdr (25) 

In general, the function g (r) depends on characteristics of the medium, 
and it can be related 

g (r) = a 2 G (e x r) , (26) 
to the so-called structure function of the medium 

G ( r ) = (E^ 5(Xi-x *+"- r) )' (27) 

Xj are coordinates of the i-th scattering center. In Eq. (126]) a 2 is the scattering 
area of the particle and e x is unit vector along the trajectory of the particle. 
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This relation can be established rewriting Eq. ( 12 7p in the form 

G(r) = J>;(r), U' n {v)dv = l (28) 

<( r ) = (5(xi-x. J+rt -r)) 

Although the sum in Eq. (|28p is going over all scattering centers in the 
medium, only centers along the trajectory of the particle enter into Eq. (126]) . 
Expanding Eq. (12~T|) in powers of ip (q) and taking the inverse Laplace trans- 
form of each term of the obtained series, we get 

g(r) = J2^n(x), fip n {r)dr = l (29) 

where ip n (r) = (5 (xi — Xi +n — r)) is the probability distribution of the dis- 
tance r = Xi +n — Xi between n consequent collisions along the trajectory of 
the propagating particle (see Eq. (fTTj) ). Comparing Eq. (T2B1) and f[2"9l term 
by term, we reproduce relation (126]) between functions g (r) and G (r). 



3.3 Continuous limit 

In the case of scattering medium with fractal dimension D the Laplace trans- 
form of the function g (r) in the long wavelenth limit aq 1 has the form: 

g(q) = (aqy a } a = D — 2 (30) 

In Appendix Owe derive asymptotic solutions of Eq. (123 p close to the surface 
and far from it. Close enough to the surface the distribution of energies will 
have large peak at A = 0, describing non-scattered particle 

f(A,x)~6(A)*(x), (31) 

where the function \l/ (x) is the survival probability that the particle does 
not scatter at the depths smaller than x. Using large p asymptotes of the 
function u (p) ~ 1/a defined in Eq. we find from Eq. ( 15T|) of Appendix O 

* Or) = E a {{x/af] (32) 

where the function E a is defined in Eq. (|46p of Appendix [Bj 
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The process with survival probability ^ (r) = E a (r) is known as the 
Mittag-Leffler renewal process, that is described by fractional differential 
equation 

d a ^ (r) 

^-ZlU. = -* r (33) 
dr a v ' v ' 

where d a /dx a is the Caputo fractional derivative 



dr a Til -a] 



f r W (x) 

/ -^dx, < a < 1. (34) 

Jo (r-x) 



It was introduced by Caputo in the later 1960s for modeling the energy 
dissipation in the rheology of the Earth [12]. 

Far enough from the surface the loss of the energy is determined by mul- 
tiple scattering processes, when we can expand 



Ai f°° 
us (jp) ~ p — , Ai = / desw 
a Jo 



Calculating the inverse Laplace transform of Eq. fl52|) in Appendix [Cj we get 



A fa 
Ai \x 



(35) 



where the Wrize function W a is defined in Eq. (H0|) of Appendix [B] 

The function / (A, x) (J35ll is the solution of the space-fractional drift 
equation of the order a: 

dfjAx) a a d a f(A,x) 

dA Ai dx a 1 ' 

This distribution function (135]) describes pure renewal process with anoma- 
lous exponent a. [9] This process can be modelled as the series of jumps of 
the energy A with the amplitude A x each happened at renewal points x sep- 
arated by random discrete intervals I distributed according to the power law: 
ip (I) ~ see Fig. ffl 

In the continuous limit we can generalize the simple model of section 13.11 
using the thinning procedure (13]. In this procedure for each positive n a 
decision is made: the scattering event is maintained with probability p or 
it is deleted with probability 1 — p. In the limit p — > the amplitude of 
scattering become smaller and smaller, their number in a given span of space 
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larger and larger, and the ballistic trajectories between scattering events 
smaller and smaller. In this limit there are no ballistic trajectories anymore 
and we come to continuous medium. 

When the energy distribution w (e) decays quicker than e~ 2 energy losses 
for the thinning procedure model are described by asymptotically universal 
Mittag-Lemer distribution. This distribution is characterized by the spatial 
scale b ^> a at which the particle is scattering with the probability about 
1/2. In the random medium with power- low correlations the parameter b 
determines the amplitude of the Laplace transform of the function g (r): 

~g(q) = (bqy\ ag«l (37) 

The Mittag-Lemer asymptotic distribution corresponds to the model of 
constant energy losses A = A\n (section I3.ip with asymptotic function 
9 (<?) (I37|) . Calculating the inverse Laplace transform of expression fl53|) in 
Appendix O and renormalizing the spacial scale a — > b we get: 



/. W = S^_£f )[(£)■]. (38) 



(xjVf 
n\ 

where is n-th derivative of the function E a , Eq. ( I46p of Appendix [Bj and 
the average over this distribution n (x) = (x/b) a /Y (1 + a) f|T7|) determines 
average losses at distance x from the surface, (A (x)) = A\n (x). At n = 
the Mittag-Lemer distribution turns to the survival probability, Eq. fl3T|) . 
while in the limit x ^> b it turns to the renormalized function (1331) . In the 
case a = 1 the distribution (1311 takes the well known Poisson form 



in (x) = — |-e" n , n = - 
nl b 

and may be considered as generalization of Poisson distribution of scattering 
events for the case of fractal medium. 

In the case of slowly decaying energy distribution at scattering events 
w (e) ~ e' 1 ^^ (0 < (3 < 1) the scattering on fractal structures has the 
form of convolution of Levy flight processes characterized by the exponent 
(3 with anomalous scaling processes characterized by the exponent a. Such 
combined process is described by fractional space-energy differential equa- 
tion, and average energy losses grow with the distance x from the surface as 
(A (x)) ~ x a IP . The effective exponent a/{3 of this process can be smaller or 
larger than 1, depending on relation between exponents a and ft. 
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4 Conclusion 



We study the loss of energy of the particle moving in fractal media and 
in the system with dynamic heterogeneities formed at the critical point of 
phase transition. We show that when the distribution of energy loss during 
collisions w (e) quickly decays with the energy e the distribution function of 
particle energies is universal and depends only on fractal dimension D of the 
medium. In the case D = 3 spacial heterogeneities only weakly affect the 
scattering process, which can be described by the classical theory. [10] In the 
case 2 < D < 3 spacial heterogeneities change the character of the scattering, 
which can be described by fractional differential equations of order a = D — 2. 
Nonlocal character of fractional derivatives (see Eq. fl3~4"j) ) reflects power-law 
correlations existing in the fractal system. We show that the loss of the 
energy in fractal medium can be described by the Mittag-Lemer renewal 
process of order a, which is fractional generalization of the Poisson process 
corresponding to the case a = 1 of absence of such correlations. 

One of the most important applications of this theory is propagation of 
particle through the percolation scattering structure with fractal dimension 
D ~ 2.49, when the exponent a ~ 0.49. Similar exponent a = 1/2 is 
obtained for fractals with lattice animals structure - a set of randomly con- 
nected sites on a lattice. [H] The lower boundary D = 2 (corresponding to 
random walk structures) of applicability of our consideration equal to the 
dimension d = 2 of the interface. Therefore in the case of fractal dimension 
D < 2 of scattering clusters the losses of the particle energy per unit area 
decrease with the rise of the interface area. 

We also derived general kinetic equation (I2"3"j) for the average distribution 
function of energy losses in random fractal medium, that can be considered as 
generalization of the Landau equation ([1]) for homogeneous medium with con- 
stant density of scattering centers, g (r) = 1/a. In heterogeneous medium the 
integral kernel g (r) of this equation is proportional to the structure function 
G (r) of the medium. The non-local character of the kinetic equation (T23]) is 
related to the presence of strong non-local correlations in a fractal medium. 
The consideration of scattering of paricles in turbulent medium needs addi- 
tional study because of multifractal structure of turbulent flows. 
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A Estimation of exponent a for fractal medium 

Here we present simple scaling estimation of the survival probability \I/ (r), 
Eq. (Q, in the case of scattering of particle in random medium with fractal 
dimension D. We consider all particle trajectories colliding with the same 
scattering center as identical. Average number of such trajectories starting 
from one scattering center can be found by draw the sphere of radius r 
around this center, see Fig. [2j The total number of scattering centers inside 
this sphere is N (r) ~ (r/a) D . Note, that scaling consideration can be applied 
to fractal objects with some care: to get the right scaling on the scale r we 
have to place all these scattering centers randomly at a distance ~ r from 
the center of the sphere. The propagating particle collides (first time, as in 
the case of the first passage problem, see Fig. [2]) only with small survival 
part \& (r) of these centers. Therefore, the total number of different ballistic 
trajectories is N (r) \1/ (r) and their total area on the sphere is a 2 N (r) \P (r). 
On the other hand, projections of these trajectories on the sphere cover the 
whole area ~ r 2 of this sphere: 

a 2 N (r) (r) ~ r 2 

Solving this equation with respect to the survival probability we find 

\]> ( r ) ~ ~ _ 

1 J a 2 N{r) \rJ 



(39) 



B Wrize type functions 

The Wrize function is defined as 

^w = E ar (i-l-ao < 40) 

We can present this function in integral form using corresponding presenta- 
tion of Gamma function for non-integer —z 



1 z 
fjz) ~ 2^ 



-ty z e- l dt (41) 



c 
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Figure 2: Scaling estimation of survival probability \I/ (r). The number of 
scattering centers inside the sphere of radius r is N(r) ~ (r/a) . The 
particle collides first time with only N (r) \l/ (r) centers. The total surface 
area of these centers a 2 N (r) \1/ (r) ~ r 2 equal to the area of the sphere. 



The contour C of integration encircles positive axis of complex variable t. 
Substituting this expression into Eq. (jlDj) and calculating the sum over I we 



Calculating the integral fj42|) by the steepest descent method, we find its 
asymptotic behavior at large z ^> 1: 



find 




(42) 



The function W a (z) is normalized by the condition 




(43) 




(44) 



e z for a = 



(*) = < ^e" 274 for a = 1/2 



(45) 



5 (z — 1) for a = 1 



The one-parameter Mittag-Leffler function is defined by the series 



oo 




E a (z) = 



r(l + la)' 



(46) 
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E a (z) = ^- [ ( l ] , a e- l dt (47) 



Its integral presentation can be found similar to Eq. 

_ r 1 

2tt J c z + (-t) 
and at large z ^> 1 and a < 1 it decays as: 

B -^r(r^ (48) 

We show also known explicit expressions for this function: 

^1 + z)~ x for a = 
E a (z) = { e z2 erfc (z) for a = 1 /2 (49) 
e _z for a = 1 

C Solution of CTRW 

The solution of Eq. ([22]) with function g (g) ( 130|) has the form: 

/(P.g)= 7 TV/ a- Q (5°) 
9 + w (p) (ag) 

where a; (p) = [1 — w (p)] /a. Calculating the inverse Laplace transform over 

q 



POO 

f(p,x)= / dAe- pA f(A,x) 
Jo 



we find in the limits of small and large x: 

a) At small x one can expand Eq. (!50l) in powers of small q~ a , and we get 
simple analytical form of the Laplace transform of the distribution function 

f(p,x) = E a [u(p)a 1 - a x a ] (51) 

where E a is the one-parameter Mittag-Leffler function, defined in Eq. (1471) 
of Appendix [Bl 

b) At large x one can expand Eq. (1501) in powers of small q a , and we find 
/(P.*) = Er(l-aO [aw (p)(x/arr ^ 



=i 
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In the case of constant losses of particle energy A = Ajn the function 
ui (p) = [l — e~ pAl ] /a, and we find from Eq. (|5TT|) 



(ag) (1 ~ a)n 



/n(g) = T ^ U - a1 n + l (53) 
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